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ABSTRACT 



Aims. We report the discovery of periodic flares of 6.7 GHz methanol maser in the young massive stellar object G22. 357+0.066. 
Methods. The target was monitored in the methanol maser line over 20 months with the Torun 32 m telescope. The emission was also 
mapped at two epochs using the EVN. 

Results. The 6.7 GHz methanol maser shows periodic variations with a period of 179 days. The periodic behavior is stable for the last 
three densely sampled cycles and has even been stable over ~ 12 years, as the archival data suggest. The maser structure mapped with 
the EVN remains unchanged at two epochs just at the putative flare maxima separated by two years. The time delays of up to ~ 16 days 
seen between maser features are combined with the map of spots to construct the 3-dimensional structure of the maser region. The 
emission originating in a single ~100 AU layer can be modulated by periodic changes in the infrared pumping radiation or in the 
free-free background emission from an Hll region. 

Key words, stars: formation - masers - ISM: individual object (G22.356+0.066) - radio lines: general 



1. Introduction 

Variability is a definite characteristic of methanol masers asso- 
ciated with massive young stellar objects (MYSO). The obser- 
vations of 6.7 GHz maser flux in several sources have shown 
a significant level of variability on timescales of a few days to 
several years ( Caswell et a l l 1995 1 MacLeo d & G aylard I 19961 
Goedhart et al. l2004l |2009l Sugiyama et aL l2008l l. The possi- 
ble causes of variability are changes in the pump rate or the 
maser path length, such as those caused by large scale motions 
(Caswell et al. |19951 l. A long-term monitoring of 54 sources has 
revealed a diversity of types of behavior at 6.7 GHz, including 
periodic flares with periods of 132-520 days found in six ob- 
jects (Goedhart et al.|2004Ji. For the best-studied periodic source 
G9.62-hO.20E, van der Walt et al. ([2009 ) propose a colliding 
wind binary that modulates the background radio continuum 
and/or pumping radiation with a period of 244 days. Araya et 
al. (2010) report the detection of ~237day periodicity in the 
4.8 GHz formaldehyde and 6.7 GHz methanol maser lines to- 
wards MYSO candidate G37.55-hO.20. As an alternative sce- 
nario for those variations, they propose regular changes in the in- 
frared pumping radiation due to periodic accretion onto a young 
binary system. 

In this Letter we report the discovery of periodic varia- 
tions in the 6.7 GHz maser flux in G22.357-h0.066. The posi- 
tion of 6.7 GHz maser (Schutte et al. l 19931 ) measured with VLBI 
(Bartkiewicz et al.|2009| coincides within 0'.'05 with the water 
maser (Bartkiewicz et al.|201 1). The source of 1.2 mm dust con- 
tinuum (Beuther et al. 2002) is offset by 2'.'9 from the methanol 
maser The extended (-23x10") HII region (White et al. l2005l ) 
and HllOff recombination line and 4.8 GHz formaldehyde ab- 
sorption (Sewilo et aI. I2004l ) are offset by more than 10", cor- 



* Figures 3 and 4 are only available in electronic form at 
|http://www.aanda.org| 



responding to 0.2 pc for the near kinematic distance of 4.86 kpc 
(Reid et al. l2009l ). Using the Torun 32 m dish data we found peri- 
odic maser flares with a high degree of stability and were able to 
determine time delays between the individual spectral features. 
Those delays and high angular resolution maps obtained with 
the EVN are used to recover the three-dimensional structure of 
the maser emission. There is convincing evidence for radiative 
excitation of maser spots by a variable central object. 



2. Observations 

The 6668.519MHz spectra were obtained with the Torun 32 m 
radio telescope from 2009 June to 2011 February at irregular 
intervals of one to six weeks. From 2010 September, when it be- 
came clear that the source was varying on shorter timescales, the 
monitoring interval was shortened to 1 -20 days. The half power 
beam width of the antenna was 5'.5, and the system tempera- 
ture was typically around 40 K. A frequency-switching mode 
was used. The signal of dual circular polarization was ana- 
lyzed in a 4 MHz bandwidth divided into 4096 channels, yield- 
ing a velocity resolution of 0.05 km s ' after Hanning smooth- 
ing. A typical rms noise level in the spectra after averaging 
the two polarizations was 0.2-0.3 Jy. The flux density scale 
was established by observations of 3C123 assuming flux den- 
sity from Ott et al. ( 119941 ), while the stability of the system was 
regularly checked by observations of the non-variable source 
G32.745-0.076 (Caswefl et al. 1995). Careful inspection has 
shown that some features in this source are essentially constant 
within 2-3% over the two years of monitoring. The accuracy of 
the absolute flux calibration was usually better than 10%. 

Archival spectra were retrieved from the Torun 32 m 
archive for MJD 1392, 2146, 4657, and 4849, where 
MJD=JD-2450000. The spectral resolution was 0.05 km s ', 
while the resulting accuracy of absolute flux was about 15%. For 
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Table 1. Variability parameters of the maser features in 
G22.357+0.066. 
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Fig. 1. Light curves of selected 6.7 GHz methanol maser features 
in G22. 357+0. 066. The 5 channel averages, labeled by the LSR 
velocity of the central channel, are shown. 



the first two epochs, the velocity scale was subject to +0.4km s~' 
error, and we assumed that the peak velocity was SO.lkms"'. 
The target was observed with the EVN at two epochs: 2007 
June 13 and 2009 May 29. Details of the observational setup, 
data reduction, and results for the first epoch are described in 
Bartkiewicz et al. (I2009I I. Uniformly weighted maps were ob- 
tained at both epochs with the respective beam 8.7x4.8 mas and 
6.0x4.4mas and spectral resolution of 0.18 and 0.09kms 

3. Results 

The light curves of the three strongest features in the spec- 
trum are displayed in Fig. [1] Most of the spectral features show 
cyclic flares usually followed by the return of the flux density 
to basically the same quiescent level (within ~0.2Jy) between 
flares. The feature at 79.5 km s"' exhibited a linear decrease 
(1.29 Jy yr"') of the quiescent flux during the monitoring period. 

We quantified the variability of maser flux using the variabil- 
ity indices Yl\ and Yl^. The first index is defined (Goedhart et 
al. l2004l) as 



y/i = ( £](m(/,) -m)^- Yjn{U) - W)^)lm 



(1) 



where is the total number of observations, m(f,) the observed 
flux density at a given spectral channel at f, epoch, n(f,) the flux 
density at any emission-free spectral channel at f, epoch, and m 
and n are the average flux densities. The second variability index 
is Vh = (5,™x-'5inin)/(5max+'Smm), whcrc 5max and 5 min are the 
highest and lowest flux densities, respectively. The spectrum of 
VIi is superimposed on the 6.7 GHz flux density spectrum (Fig. 
12^). The values of VI\ and VI2 for the most recognizable features 
are given in Table [T] We notice that the blue-shifted emission, 
i.e. at velocities lower than the systemic velocity of 84.2 km s"' 
(Szymczak et al. l2007l l. is significantly more variable than the 
red-shifted emission. 

Figure |2j; shows the light curve of the 80.1kms"' feature. 
The new data are supplemented by the archival 32m telescope 
observations at four epochs (see Sect. 2). No velocity shift of the 
strongest feature (80.1 kms"') was detected within the velocity 
resolution for ours and the Walsh et al. (.19971 data. The flux 
variations are fitted using the equation (David et al. 11996) : S (t) - 
j^s(t) _|. where A and C are constants and 5(f) - [{bcos(a>t + 
(f>))/(l - /sin(tijf + ip)] + a. Here, b is the amplitude measured 
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Fig. 2. a) Average 6.7 GHz maser fine profile of G22.357+0.066 
(blue line) with the variability index VIi superimposed (red dot- 
ted line), b) Folded flare profile at 80.1 kms ' feature. The dot- 
ted line represents the best fit to the curve. The symbols mark 
different cycles, c) Time series for the 80.1 kms"' feature with 
the best-fit curve superimposed. The arrows mark the dates of 
the EVN observations. The residuals from the fittings in b) and 
c) are shown by the crosses. 



relative to the average value a, oj - 2n/P, P is the period, the 
phase, and / the asymmetry parameter defined as the rise time 
from the minimum to the maximum flux, divided by the period. 
The best fit to the observational points is shown in Fig|2};, and the 
inferred period is 179.2+0.6 days. The time of flux maximum of 
the last cycle was MJD=5518.8±1.6. Smafl residuals (<2.2Jy) 
for the 1999 and 2000 points suggest the high regularity of flares 
on a timescale of ~12years. Furthermore, a peak flux of 13 Jy 
reported for 1994 April 8-12 observation (Walsh et al. 1997 1, i.e., 
16-20 days before a back-extrapolated maximum supports this 
suggestion. 

Figure|2j) shows the shape of the flare at feature 80. 1 km s"' . 
The data are folded on the precise 179.2-day period, and the ori- 
gin of time axis is arbitrary. The rise and decay times of the flare 
were 14.5 + 1.7 and 42.9+1.7 days, respectively. The full width 
of flare at half maximum (FWHM) was 19.6±0.7days. The val- 
ues of FWHM for other features are given in Table [1] There is 
a slight increase in the flare duration for features with higher 
radial velocities. This tendency is also visible in the data from 
the last cycle (Fig.O; for the blue-shifted emission; lower than 
80.5 kms"', the flare lasted on average 65.8+2.6 days and in- 
creased up to ~80days for 81.01 and 8 1.45 kms"' features. 

The most striking characteristic of the flare are the time de- 
lays of the maximum flux between features of different velocities 
(Fig. O. To measure the time delay precisely, we used the dis- 
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Crete autocorrelation function (Edelson & Krolik [T988] l. because 
the light curve is covered with unevenly sampled data. The cal- 
culated time delays. At, for the central velocities of maser com- 
ponents are listed in Table[T] A negative lag means that the maser 
feature flares before the reference feature (80.1 kms"'). The de- 
lays range from -3.4 to 13.3 days. The estimated uncertainties in 
the time delays are usually below 0.7 days but are higher than 
1 .5 days for the components, marked in italics, that have a low 
signal-to-noise ratio and small flare amplitudes. 

The 6.7 GHz maser structure of the target observed with the 
EVN at MJD 4981, i.e. exactly at the time of the putative max- 
imum of the flare (Fig.|2]3), is shown in Fig.|4] For comparison, 
we present the map obtained by Bartkiewicz et al. (2009 1. Their 
observations were also carried out just at a flare maximum (MJD 
4264), as estimated from the 32m dish data. In general the struc- 
ture of the maser emission is well preserved over two years, as 
the relative positions of the spots are the same within ~0.5 mas. 
There are only two spots for which the morphology has appar- 
ently changed. The cluster of spots centered at 84.86 km s"' was 
more scattered at the first epoch. No emission from the S-E part 
of this cluster was detected at the second epoch. Similarly, the 
emission from the northernmost part of the cluster of spots near 
88.47 km s"' disappeared after two years. These two red-shifted 
clusters show relatively weak emission (<1.7 Jy). Their flare pro- 
files are flatter and broader than those of the blue-shifted emis- 
sion (Fig. [3]). One can suppose that they are diffuse and have low 
brightness, which can be filtered out more at the second epoch 
when observed with a beam 1.6 times smaller and with spec- 
tral resolution 2 times higher than at the first epoch. No new 
maser spots are formed over two years. The brightness of indi- 
vidual maser spots of the blue-shifted emission usually changes 
by less than 1 -2 Jy, which is within the error range of absolute 
amplitude calibration. We conclude that there are no significant 
differences in the source structure at the two epochs. 

4. Discussion 

The detection of periodic variations in the 6.7 GHz methanol 
maser flux in G22. 357+0. 066 extends a class of periodic maser 
sources (Goedhart et al. l2004l l. We note that the flare profile of 
the target is very similar to what is reported in G9.62+0.20E 
at 12.2 GHz (van der Walt et al. |2009 |. Specifically, the ratio of 
the rise to the decay time is 0.34 in G22. 357+0.066 and 0.37 
in G9.62+0.20E, and the respective values of the relative ampli- 
tude of flux density variation are 1.5 and 2.0. This suggests that 
the same underlying processes drive the periodicity in the two 
sources. 

The cause of the periodic flares of methanol maser sources 
is still unknown. Several possibilities have been considered 
(Goedhart et a l.l2009l van der Walt et al. l2009l Araya et al. l20T0l 
van der Walt '201 1\ including (i) disturbances of the masing 
region by shock waves or clumps, (ii) variations in the back- 
ground radio continuum and/or of pump photons due to stel- 
lar pulsations, modulations of the radiation from a colliding 
wind binary, or periodic accretion of circumbinary disk mate- 
rial. For G22. 357+0. 066 we found that all the features return 
their 6.7 GHz flux to basically the same quiescent level between 
flares and that the velocities of different features and the spec- 
trum shape are stable over at least 12 years. Moreover, the spa- 
tial structure on a mas scale remains essentially the same for the 
two flares spanning two years. These facts imply that the velocity 
and spatial structure of the source remain unaffected by whatever 
mechanism underlies the periodic variation. Furthermore, the 
maximum projected separation of 1050 AU between the maser 



spots in the source, for the distance of 4.86 kpc, is well within a 
range of typical distances between the maser region and the driv- 
ing MYSO (Bartkiewicz et al. 2009). If the flares are triggered by 
a propagating density wave or shock front, its velocity should be 
higher than lO'* km s"' to account for 179 day periodicity. Such a 
fast shock would easily dissociate the molecules (Hollenbach & 
McKee 1989). Therefore we exclude shock waves as a possible 
mechanism of the flaring. 

The time delays between flares of individual features, with 
the highest values of ~16 days between the flare in the blue- 
shifted feature at 79.5 km s"' and the extreme red-shifted fea- 
ture at 88.5 km s"', clearly implies the radiative coupling of a 
single triggering process with the maser regions. This is con- 
sistent with very similar flare profiles for the emission features 
which arise from single cluster of spots, i.e. at velocity lower 
than ~81kms However, the flare profile of the 81.5kms ' 
feature is different (Figs.[Tl[3]l. This feature is a blend of emis- 
sion from the two spatially separated clusters (Fig. |4), and its 
blurred flare profile is likely due to the different time delays of 
each cluster The red-shifted features near 85.0 and 88.5 km s ' 
have flat and poorly pronounced flare profiles. Assuming that the 
maser flaring is caused by changes in the infrared pumping radi- 
ation, the red-shifted emission may come from the far side of the 
envelope and it is less bright and less saturated due to the lack 
of background radiation. Alternatively, it can be due to either 
strongly directional amplification of the maser emission, where 
backward propagation dominates, or blending of low intensity 
extended emission. However, if the flaring of all the maser fea- 
tures is due to changes in the free-free emission of a background 
HII region, then all the maser features must be on the near side 
of the HII region and no simple geometry of the maser region 
can be recovered from the time delay measurements. 

We used the time delays and the EVN map to construct the 
three-dimensional structure of the maser region. For a scenario 
of modulation of the infrared pumping, we assumed that a source 
triggering the periodic variations coincides with the line of sight 
to the strongest feature at 80.1 km s"' and that the maser emis- 
sion arises in a spherical envelope of mean radius ro. The value 
of ro is calculated by minimizing residuals of the least-square fits 
of the envelope to the positions of spots. The best fitted value of 
ro is 710+140 AU, and the maser spots lie in a layer of ~100 AU 
thickness inclined 47° from the line of sight (Fig.|5]l. The over- 
all size of the maser structure is ~1 500x600 AU. It is a some- 
what surprising result, which may indicate that the maser emis- 
sion forms in a circumstellar disc or torus. We note that this is 
not unique model because neither the location of the triggering 
source nor details on the pump rate variations are known. 

The strongest maser spot coincides within 0'.'62 with the 
GLIMPSE source G022.3564+00.0662 (Fazio et al. l2004] i: i.e., 
the total emission seen by the EVN falls within one Spitzer 
pixel of nominal size ~l'.'2. The colors of this counterpart source 
are typical of an MYSO (Cyganowski et al. l2009l l. No compact 
5 GHz radio continuum (lcr=0.34mJy b ') or any H2CO ab- 
sorption was detected with the beam 1'.'6 x I'.'l (Araya 2010, 
priv. comm.). The extended (~23"xl0") HII region (White et 
al. 2005) and HI lOo- recombination line and 4.8 GHz formalde- 
hyde absorption (Sewilo et al. l2004l l are offset by more than 10", 
while a faint (1.02mJybeam ') and extended 8.4 GHz emission 
is located by 14'.'8 to the N-E of the source (van der Walt l2003) . 
corresponding to 0.2-0.3 pc. There is a suggestion that the maser 
lies at the edge of an evolved HII region of ~5' size that may 
trigger a new generation of stars (van der Walt et al. 2003]- 
Therefore it is probable that more sensitive observations will 
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Fig. 5. Three-dimensional view of 6.7 GHz maser in 
G22. 357+0. 066. The z axis is towards the observer The 
symbol sizes are proportional to the variability index VI2 (Table 
[T]), and the colors mark the LSR velocity scale as in Fig. |4] The 
bars mark the distance of spots to the plane shown by the grid. 



identify a source triggering periodic variation of the methanol 
maser 

A model of a colliding wind binary was proposed to explain 
the best-studied periodic source G9.62+0.20E (van der Walt et 
al. 2009; van der Walt 2011). It provides the mechanism of pe- 
riodic changes in the background radio continuum and infrared 
radiation. Regularly modulated pulses of ionizing radiation pass- 
ing through a volume of partially ionized gas against which the 
maser is projected produce periodic changes in the electron den- 
sity so that the background free-free emission is variable. Within 
the frame of that model, the observed decay time of the 12.2 GHz 
maser flare in G9.62+0.20E is consistent with the characteris- 
tic recombination time of a hydrogen plasma with densities of 
10^ - 10* cm""*. As there is a hint of an extended radio con- 
tinuum source in the field in the case of G22. 357+0. 066, it is 
reasonable to apply the van der Walt at al. model. We fitted eq. 3 
(van der Walt et al. |2009l l to the decay part of the 6.7 GHz flare in 
the source (see Fig. 2b). The best fit was obtained for the equilib- 
rium electron density ranging from5xl04 to 1x10^ cm and the 
electron density at the decay start time of 1.03x10* cm"-'. These 
values are well within the ranges reported for G9.62+0.20E (van 
der Walt et al. l2009l l where two components of the decay of the 
flare profile were fitted. This suggests that the maser flare in 
G22. 357+0. 066 can be due to changes in the number of free- 
free background photons from an HII region. 

Araya et al. (2010^ propose that periodic heating of the dust 
and increasing the infrared radiation can be caused by accretion 
of material from the circumbinary disk onto the protostars or 
accretion disks. As pointed out in van der Walt (.201 1) . this sce- 
nario is less probable for periodic masers, because for a system 
with an eccentricity of 0.1, a significant time lag is expected in 
the accretion rate onto the two stars, and variations in the to- 
tal accretion rate will not show the sinusoidal-like pattern. For a 
highly eccentric system, the accretion rate changes rapidly near 
periastron (Artymowicz & Lubow [T996) . Although the scenarios 
of the colliding wind binary and periodic accretion onto a young 
binary system provide a periodic source of photons heating the 
circumstellar dust, thereby aff'ecting the pumping infrared radi- 
ation field, the typical cooling time for the optically thick case 



of a few days (van der Walt et al. l2009l l is too short to explain 
the decaying part of the maser flare profile. Monitoring of the 
source in the radio and infrared will provide major progress in 
explaining the causes of its periodic variability. 
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Fig. 3. 6.7 GHz flux density of G22. 357+0.066 as a function of time and velocity. The timescale covers the last cycle of variability 
from September 2010 to February 201 1 . 
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Fig. 4. Maps of 6.7 GHz metha nol ma ser in G22. 357+0.066 obtained with the EVN at two epochs; MJD 498 1 (a) and MJD 4264 (b, 
adapted from Bartkiewicz et al . 120091 1. The coordinates are relative to the brightest spot. The sizes of symbols are proportional to the 
logarithm of the brightness of maser spots. The colors of circles relate to the LSR velocity scale presented in the wedge. The groups 
of spots are also labeled by their central velocities or the velocity ranges. The linear scale, derived for the distance of 4.86 kpc, is 
shown by the horizontal bar. 



